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We investigate pairwise correlation properties of the ground state (GS) of finite antiferromagnetic 
(AFM) spin chains described by the Heisenberg model. The exchange coupling is restricted to nearest 
neighbor spins, and is constant Jo except for a pair of neighboring sites, where the coupling Ji may 
vary. We identify a rich variety of possible behaviors for different measures of pairwise (quantum and 
classical) correlations and entanglement in the GS of such spin chain. Varying a single coupling affects 
the degree of correlation between all spin pairs, indicating possible control over such correlations 
by tuning Ji. We also show that a class of two spin states constitutes exact spin realizations of 
Werner states (WS). Apart from the basic and theoretical aspects, this opens concrete alternatives 
for experimentally probing non-classical correlations in condensed matter systems, as well as for 
experimental realizations of a WS via a single tunable exchange coupling in a AFM chain. 
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I. INTRODUCTION 

Quantum behavior at the microscopic scale is well es- 
tablished both theoretically and experimentally. Strictly 
quantum phenomena observed so far have confirmed the 
validity of quantum mechanics postulates, in particu- 
lar its non-local character^ Entanglement is certainly 
one of the most intriguing and basic of such quantum 
phenomena, as it is a fundamental resource in perform- 
ing most tasks in quantum information and quantum 
computation^ Recently it has been realized that entan- 
glement is not the only form of quantum correlations, 
and that separable quantum states, that were expected 
to show strictly classical behavior, in fact can exhibit 
quantum correlations other than entanglement. A mea- 
sure of quantum correlations that has recently received 
much attention is the quantum discord^ It has been an- 
alyzed both theoretically^— and experimentally*^— 

Experiments demonstrating and quantifying entan- 
glement have mostly dealt with photons^ — and 
atoms i 17 : 18 Solid state systems could in principle pro- 
vide more compact and stable hosts to study such ef- 
fects. It is also desirable to demonstrate control over 
pair-related quantum phenomena in solid state systems, 
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or 



such as spins of electrons bound to quantum dot: 
to impurities^ in semiconductors. 

We investigate here pairwise correlations and entangle- 
ment in small (few sites) magnetic chains. The magnetic 
moments are the spins of electrons individually bound 
to a short linear array of quantum dots in GaAs or Si- 
based nanostructures, or of donors in Si. Neighboring 
spins interact via exchange coupling. An important as- 
pect of such fabricated chains is the tunability^ of the 
coupling between spins in specified pairsi^^r— Changes in 
this coupling affect the pair correlations, and we identify 
general trends to be expected when one coupling varies in 
such chains. As a general trend obtained, the behavior of 
the different correlations studied is similar, seeming to il- 
lustrate a monogamous nature^— , even for correlations 



other than entanglement. Also, we show that a class of 
two-spin states of even-numbered chains are exact real- 
izations of Werner states (WS). These states are useful as 
they simulate the effect of the environment in destroying 
coherence and entanglement in a quantum system. 

This paper is organized as follows. Our model is pre- 
sented in Sec. [H] where the Hamiltonian and the den- 
sity matrix are introduced. Measures quantifying corre- 
lations and entanglement used here are briefly reviewed 
in Sec. IIII1 In Sec. HVI the behavior of small model chains 
is presented. In Sec. [V] we obtain and discuss the realiza- 
tion of WS in our spin model. Our final conclusions are 
given in Sec. [VD 



II. MODEL HAMILTONIAN AND TWO-SPIN 
DENSITY MATRICES 

We consider the isotropic Heisenberg model for anti- 
ferromagnetic (AF) spin-half finite chains. The Hamilto- 
nian for this system is given by 



A' 



H — ^ ^ Jij S{ • Sj , 



(1) 



■;3) 



where the summation is over first neighbors of an open 
chain with N spins, and the exchange coupling is positive 
(Jij > 0). We obtain the spectrum for the Hamiltonian 
in Eq. (JT]) in the product basis of the Sf eigenstates, of 
the form \sis 2 ■ ■ ■ sjv) = |si) ® |s2) (8) • • • ® |sjv) where 
Si represents up or down spin at site i, i.e., |s,) = {|t 
)i, \l)i}. The total spin operator in the z-dircction, S J = 

X)i=i Sf s i s a good quantum number, as [H, Sj] = 0, 
and so the Hamiltonian exhibits a block-structure, each 
block having a well-defined value of S'J. For AF coupling, 
the lowest state is in the lowest Sj subspace, so even-N 
chains have a non-degenerate S J = ground state (GS). 
Odd-A^ chains GS is doubly-degenerate, one state from 
each block Sj = 1/2 or = -1/2. 
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The simplest AF system is that of two spins, for which 
the GS is a singlet: 



\s) = ^= (in) - nt» 



(2) 



This is a maximally entangled state, in the sense that 
knowledge of the global state of the two spins gives no 
information about the state of either spin, while by choos- 
ing an observable, measurement of the state of either spin 
gives with certainty the state of the other. This quan- 
tum correlation is stronger than could ever be possible 
classically. 

The two-spin density matrices pij for spins located at 
sites i and j, are obtained from the full density matrix 
Po = \ipo)(' l Po\ by tracing out the degrees of freedom of 
the other TV — 2 spins. These are in general mixed states, 
written as a weighted sum of pure states: py = ^2iPipi, 
where pi is the weight of pi = \ipe)(ipi\. 

It has been shown that the symmetries of the Hamilto- 
nian (JlJ imply that the two-spin density matrix for any 
pair of spins in the S z basis {|ft)i |t4)) lit), has the 

general form^^ 



( a 

h z 

z b 2 

V d 



(3) 



with Trpij = 1, &1&2 > \z\ 2 and ad > 0.— The single spin 
density matrix pi corresponding to spin i is obtained by 
tracing out spin j. 



III. QUANTIFYING CORRELATIONS 

Several quantities describe correlations in quantum 
systems; we summarize below those adopted in our study. 
The amount of information contained in a quantum sys- 
tem, assumed here to be a spin array described by p, 
is traditionally measured through the Von Neumann en- 
tropy, 



S(P) = -Trplog 2 p , 



(4) 



which is the quantum generalization of the classical Shan- 
non entropy. When spins i and j interact, they usually 
share some information, and this shared information may 
be quantified by the quantum mutual information (QMI), 



l(pij) = S(pi) + S{pj) - S(pij) , 



(5) 



which includes all (quantum and classical) correlations 
between spins i and j. QMI may also be defined as 



J(pij) = S(pi) -S(pi\i) 



(6) 



where S(piy) is the quantum conditional entropy, i.e., the 
entropy of spin i when the state of spin j is known. Con- 
ceptually, obtaining S(pi\j) involves measurement over 



one spin, say j. Restricting ourselves to projective mea- 
surements, denoted by the set of projectors {n^}, each 
possible state of the two-spin system after the measure- 
ment is of the form 

Pm = — (i, ® rp m ) P4J (i, ® n4) , (7) 

Pm 



where 



p m = Tr(l i ®Ili n )p ij (l i ®ISl n ) (8) 



is the probability of outcome m and lj is the identity 
operator for spin i. Quantum conditional entropy is then 
defined as a weighted sum of the entropies associated to 
the possible states p m , 



VmS{p m ) 



(9) 



Classically, T and J are equivalent, but the quantum gen- 
eralizations may differ. Their difference defines Quantum 



S(pij) =l(p ij )-J{p ij ) , 



(10) 



which also depends on the measurement basis. Note that, 
contrary to I, J depends on the measurement basis used 
to specify the state of spin j and is not symmetric with 
respect to interchange i -H- j. Quantum correlation (Q c ) 
is defined as 



f(Pij) = min &(Pij) - 
{nu 



(11) 



where the minimization is with respect to all possible 
measurement basis sets. This minimization amounts to 
finding the measurement basis that minimally disturbs 
the system, in particular extracting information about 
spin j with minimum disturbance on spin i. Defining 
classical correlation (C ) as^~— 



c °(.Pij) = max J(Pij) i 
{nu 



(12) 



it is possible to decompose the total correlation as a sum 
of quantum and classical contributions: 



1{pi 



\p l3 ) + c c ( Pij ) 



(13) 



The individual components Q c and C c are not sym- 
metric in i «-> j, since J is not, while their sum X is 
symmetric. 

Finally, the degree of entanglement between two spins 
is also investigated here through the usual measure, the 



concurrence 



ith 
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C = max{0, A} 



A 



(14) 



(15) 
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The four numbers Ai=i r .. 4 are the eigenvalues in decreas- 
ing order of the operator Pijpij, with 

Rij = ((Ty ® Oy)p i] ((Jy <g) Uy) , (16) 

where <j v is the second Pauli matrix. Concurrence has the 
property < C < 1, with C = indicating that the state 
is separable while C = 1 for a maximally entangled state, 
such as a two-spin singlet. Regarding spins interchange, 
we note that C is symmetric: C%j = Cji. 

The concurrence for systems described by the density 
matrix of Eq. (J3j) can be directly obtained from Eq. (fT4]) 
and is given by 

C = 2max(0, \z\ - Vad) . (17) 

Taking for example the singlet state of Eq. ((2|), we have 
I(= 2) = C c {= 1) + Q c (= 1) and C = 1. These values 
are extreme (maxima), meaning that the singlet is max- 
imally correlated, in the classical and in the quantum 
sense. 

In what follows, when there is no ambiguity, we refer in 
general to all quantities C, I, C c and Q as correlations. 

IV. GENERAL TRENDS 

We analyze here the behavior of correlations for two- 
spin states of even- and odd- TV linear chains as a func- 
tion of the exchange coupling parameter £ = J23/J12 = 
Ji/Jo- Explicitly we treat in this section the cases 
of N = 3 and N = 4, as depicted, respectively, in 
Figs. QJa) and [TJb) . In these chains, the coupling Ji 
varies, while the others are kept constant and equal to 
J . For N = 4 we consider also the more general case 
where J34 = J2 7^ Jo- Exact diagonalization of these 
chains were obtained here both analytically and numeri- 
cally function of £. 



{|Ttl), 14-tt)}, an d the general form of the GS state 

vector is given by 

l4 3) ) = ai|ttl> + «a|nt> + 03|ltt) • (18) 

Above, a2 is the amplitude for the classical Neel state, 
and a\ and 03 are the amplitudes of the quantum fluc- 
tuations. For J\ — Jq we obtain ol\ = a-^ = — l/v6 and 
012 = By explicit calculation it is possible to obtain 

the two-spin reduced density matrix for the state vector 
above. We have d = at Eq. [3l for two spins pointing 
down is not allowed in this subspace [for S J = —1/2, 
a = in Eq. ([3])] . The general form of the density matrix 
may be written as 

^•=PilttKttl+P2l*WI, (is) 

where 

\<f>) = Cl \n)+c2\tf) (20) 

is a general singlet, with |ci| 2 + | C2 1 2 = 1- This state has 
the property that it is maximally correlated only when 
l c i| = l c 2|- This implies that the two-spin state in the 
chain is a mixture of an entangled component with 
a separable one, |tt)(ttl: the relative weights given by 
Pi and p2- Fig. [2] shows I as a function of £ for the con- 
figuration of Fig. QJa). As Ji increases, Z12 decreases, 
while I 2 3 increases; they cross at Jj = J as, by sym- 
metry, pairs 1 — 2 and 2 — 3 become equivalent at this 
point. As regards to the second-neighbor spins 1 and 3, 
they also get correlated, although only nearest neighbors 
are exchangcd-coupled. This is an effective interaction, 
mediated by the central spin 2, as discussed at Ref. [32| 
within the context of long distance entanglement. For 
Ji <C Jo, behavior of pair 1 — 3 follows that of the pair 
2 — 3, and for J\ S> Jo it approaches that of pair 1 — 2. 
Also, I13 exhibits a maximum at J\ = Jq, corresponding 



A. Three-spins chain 

As the GS of a three-spins chain is twofold degener- 
ate, we need to choose one of the subspaces, say that of 
S J = 1/2 (considering units in which h = 1). The two 
degenerate states can be separated by an applied mag- 
netic field. The basis for the Sj = 1/2 subspace reads as 



Jq Ji 

o — 0---0 

1 2 3 

fa) 



Jo J\ Jo 

o — o — o — o 

12 3 4 

(b) 



FIG. 1. Simple spin configurations with (a) three and (b) four 
spins, with a variable coupling between spins 2 and 3. The 
other couplings are fixed to Jq. 




FIG. 2. (Color online) QMI as a function of f = J1/J0 for the 
three-spins system of Fig. [Ha). The full, dashed and dashed- 
dotted lines refer respectively to the pairs 1 — 2, 2 — 3 and 
1-3. 
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to the tendency of these two spins of forming an entan- 
gled triplet state (|tl> + ||T))/\/2, as can be verified by 
explicitly constructing pi 3 at Ji = J . The correspond- 
ing C c , Q c and C, are presented and compared to I in 
Fig.[3l where we see that the overall behavior of these four 
correlations is qualitatively the same. These results show 
that correlations tend to concentrate in the most strongly 
exchange-coupled pair, as explained by the following ar- 
gument. Two isolated spins get maximally correlated for 
any finite value of the AFM exchange coupling. When 
they are interacting with additional spins, all neighboring 
pairs will tend to form a singlet, but this is formally not 
possible. The most strongly coupled pair forms a state 
with a larger singlet character than all others, thus be- 
coming the most correlated, although not maximally^ It 



l 1 1 1 1 1 1 r 
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FIG. 3. (Color online) Total, classical and quantum corre- 
lations and concurrence as a function of £ = Ji/Jo, for the 
three-spins system at Fig. QJa). For comparison, we show 
(dotted line) the same data as in Fig. [2] for I. C, Q c and C c 
are plotted respectively in full, dash-dotted and dashed lines, 
for pairs (a) 1 - 2, (b) 2 - 3 and (c) 1-3. See text for details. 



becomes clear from these results that whenever two spins 
are maximally correlated, this precludes any correlation 
with another third spin, so that in general, as long as two 
spins share some amount of correlation, their correlation 
with a third spin is limited. This result is related to the 
so-called monogamous nature of entanglement^-"— and 
the present results seem to indicate that it applies to all 
other (quantum and classical) correlations. 

Correlations in this three-spins chain exhibit properties 
associated with the symmetry of the system which are 
easily verified by inspection of the spin arrangement and 
are also obtained algebrically. Namely pi 2 (0 = P32 
P2i(0 = P23(V£)> Pi3(0 = P3i(V0- A11 symmetries 
of p are readily transfered to all correlation functions 
described here. We recall that C and I are symmetric 
with respect to spins interchange, while p, Q c and C 
in general are not. Finally, we remark that in all cases 
we have C c < Q c < C. 



B. Four-spins chain 

We perform a systematic study of pair correlations and 
entanglement for the particular N = 4 chain in Fig.QJb), 
where J12 = J34 = Jo and the variable coupling is J23 = 
J 1. This keeps the same number of free parameters as 
the N=3 case. We discuss the more general situation 
J12 ^ J34 at the end of this section. 

Different from the N — 3 case, here we have spin inter- 
change symmetry for p (pij = pji) and for all correlation 
functions. It is possible to show by explicitly writing the 
full GS density matrix on the = four-spins basis, 
that the reduced density matrices P12 = P34, P13 = P24 
and pu = P23 for all £. As a consequence, all correlations 
follow the same symmetry properties. Regarding the 
£ l/£ correspondence given above due to the symme- 
try of the TV = 3 chains, here we find algebrically a mod- 
ified variable transformation, namely pn(£) = P23(4/£) 
which is not obvious to anticipate from symmetry argu- 
ments and, again, is transferred to all correlations. 

Fig. |4] shows correlations as functions of £ for pairs (a) 
2 - 3 and (b) 1-3. As in the N = 3 case (Figs. GUb) 
and[3^c) respectively), the overall behaviors of C, Q c and 
C c arc similar, and each pair shows C (£) < Q (£) for 
all £. However, here an entanglement/disentanglement 
transition is obtained at special values of £, while for 
N = 3 the pairs entanglement vanishes only at limiting 
situations. We have C23 = for < £ < 1, and for 
£ > 1 the pair becomes entangled. From this point on, 
C23 increases rapidly, reaches the value 1/2 at £ = 2, and 
asymptotically tends to unity. This limit corresponds 
to a singlet formed by spins 2 and 3. Independently, 
although in parallel, spins 1 and 4 also tend to form a 
singlet, as anticipated by the property /023(C) = Pia{Q- 
The other property mentioned, vis., the correspondence 
£ <H> 4/£ between pairs 1 — 2 and 2 — 3 permits one to 
infer the behavior of C12 from that of C23: C12 is unit 
at £ = 0, and decreases to 1/2 at £ = 2. It vanishes for 
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FIG. 4. (Color online) Correlations as a function of £ = Ji/Jo 
for selected spin pairs of the GS of the four-spins chain de- 
picted in Fig. [ljb) ■ ( a ) The full, dash-dotted and dashed lines 
are respectively for C, Q c and C c of pair 2 — 3. The inset 
shows the region < £ < 2. (b) Dotted, dash-dotted and 
dashed lines refer to total (J), quantum (Q c ) and classical 
(C c ) correlations for the pair 1 — 3, whose concurrence is 
zero in the full range. See text for details. 



e>4. 

As regards to the other correlations, Qf 2 (= Q^ 4 ) is 
unit for £ = 0, and decreases gradually, tending to zero 
asymptotically for large £, and Q% 3 (= Q14) increases 
from zero at £ = and asymptotically tends to one. 
In both cases Q remains non-zero in regions in which 
concurrence vanishes, indicating that separable states 
present quantum correlations not captured by concur- 
rence. A similar description holds for Cyx (= C34) and 
Cg{= Cg) [See Fig.mja)]. Here we also have C c < Q c , 
but C can be less or greater than C c and Q c . Fig. |U(b) 
shows that total, classical and quantum correlations are 
nonzero for the pairs 1 — 3 and 2 — 4 and exhibit a max- 
imum at £ = 2, while C13 = C24 = over the full range. 

In order to test the robustness and generality of our re- 
sults, we consider the case in which J34 = J2 7^ Jo- This 
is a likely experimental situation due to disorder. We find 
that there is always an entanglement/disentanglement 
transition, but at different values of £ = Ji/Jo, according 
to the value of J2. Our numerical results are consistent 
with the following transition points t;c for the pairs 1 — 2 



e c ? = 2(1 



(1 



h J34/ J12) , 

J 34 /Ji 2 )/2 . 



(21) 
(22) 



Using the values J 12 = J34 = Jo we obtain = 1 [as 
indicated in Fig. [3], and J i2 = J34 = Jo, = 4. 

We remark that symmetries discussed above, consid- 
ering J12 = J34 = Jo, remain valid in the general case 
where J34 = J2 7^ Jo- 



WERNER STATES 



As noted in Scc. lIVBl an entangled/disentangled tran- 
sition may occur for specific pairs for the N = 4 chain 
at finite values of £, e.g, C23 in Fig. Ufa) vanishes for 
< £ < 1 but not for £ > 1. As this behavior resembles 
that of the so-called Werner states (WS)piL we investi- 
gate here two-spins states in even-numbered chains to 
establish if and how they relate to WS. A WS is a one- 
parameter state that represents a mixture of a maximally 
entangled state with a totally mixed component propor- 
tional to the identity matrix.— ~— We consider here the 
following description, suited to the case at hand: 



pw 



1-p 



i 4 + p |*-)(#- 



(23) 



where I4 is the identity matrix in the four-dimensional 
space, |*_) = (|01) - |10))/V2 is one of the Bell 
states^ analogous to the singlet state of Eq. @, and 
the parameter p (0 < p < 1) gives the weight of the 
entangled component, driving an entanglement (jp > 
1/3) /disentanglement (p < 1/3) transition. 

The GS of even-TV AF chains considered here is a 
5j = state, for which the two-spins density matrix el- 
ements pij in Eq. ([3]) may be written in terms of a single 
parameter^ namely 



^ij i bi — b 2 — — — Tij , z — 2Fjj , 



where 



\S* S~. 



(24) 



(25) 



is the spin-spin correlation function. The relations (|24|) 
and (|25|) above are direct implications of rotational and 
time reversal symmetries of the Hamiltonian. Eq. ([3|) 
implies 



^./Vcvcn 



i + 4r 



1 4 -4T\S)(S\ 



(26) 



where \S) is the singlet state from Eq. ((2|). From the def- 
inition of C in Eq. {TTJ) and from the bounds — 1/4 < T < 
1 /4, it follows that when T > —1/12 the concurrence van- 
ishes and the state is separable. For T < the density 
matrix has the same structure given by Eq. (|23p . with 
the weight parameter p = — 4r, preserving the bounds 
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< p < 1. Thus, the density matrix pij of a spin pair 
of an cven-iV AF chain GS is a genuine WS when T < 0. 
Given that T is a function of the coupling ratio alone, 
r = /(£), we may conclude that, in principle, control- 
ling the coupling J\ allows control over the fundamental 
parameter that characterizes a WS. 

A. Three and four-spins chains 

As an illustration of the criteria above, we show in 
Fig. [5] the spin-spin correlation functions for the four- 
spins chain of Fig [Ijb). The correlations are always 
negative for antifcrromagnctically coupled pairs, namely 
ri 2 (= r 34 ) < and r 23 (= ri 4 ) < 0, so these pairs 
constitute representations of WS. On the other hand, 
ri 3 (= T 2 4) > 0, representing a general separable state 
(C = 0). This means that, for the S*J = states con- 
sidered, the effective coupling between second-neighbors 
spins is ferromagnetic, leading to positive T and null en- 
tanglement. This condition is not valid in general for 
odd-W chains (Sj = ±1/2 states), as for example the 
three-spins chain where although the spins 1 — 3 are cou- 
pled via a ferromagnetic effective interaction (Fi3 > 0), 
they are (weakly) entangled (Ci 3 ^ 0) for any £. 

When the relation T = /(£) is invertible for a pair 

1 — j, as is the case for r i2 and T 23 for N = 4 (See 
Fig. [5|), correlations for i — j may be explicitly obtained 
as functions of T. This leads to the result that, for pairs 
with r < 0, correlations dependence on T is universal. 
For example, Ci 2 and C 23 [Fig. Ufa)] collapse into^ 

C = 6max{0,-r- 1/12} , V < , (27) 

given by the full line of Fig.[6j The quantum and classical 
correlations for these pairs, represented in Fig. [6] are also 
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FIG. 5. (Color online) Spin -spin correlation function, T^ , 
as a function of the coupling ratio for the four-spins chain 
in FigQJb). We find Fi3 > 0, C13 = 0, in agreement with 
Fig. |4jb). The horizontal dashed line, given by F = —1/12, 
crosses Fi2 and F23, respectively, at the transition points 
£j? = 4 and £i 3 = 1. 




-r 



FIG. 6. (Color online) Correlations and concurrence as a 
function of minus the spin-spin correlation function in the 
r < region for a four-spins chain. The curves are universal 
for pairs 1 - 2, 2 - 3, 3 - 4, and 1 - 4. 

Jo J\ Jo Jo ■■ ■ Jq ■• ■ 
O O O O O O 

1 2 3 4 5 ■ • ■ N 

(a) 




FIG. 7. (Color online) (a) Schematic representation of spin 
chains considered in the paper. The "surface" spin is labeled 
"1" and all couplings between nearest neighbors equal Jo, ex- 
cept for J23 = Ji . (b) and (c) show concurrence as a function 
of the negative of the spin-spin correlation function, T, for 
odd-iV chains as represented in (a), (b) Results for the 1 — 2 
pair for N — 3, 5, 7 and 9. (c) Results for the 2 — 3 pair; the 
inset shows the rectangular dotted region. For comparison, 
the result for the WS is also shown in Figs. [7£b) and (c). 
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universal: Q^ 2 and Q 23 [Fig.[4ja)] collapse into the dash- 
dotted line of Fig. [6j while C^ 2 and into the dashed 
line. We remark that for the more general case where 
J34 = J2 =/= Jo, the two-spins states analyzed are still 
exact realizations of WS. This is a very interesting point 
in an experimental point of view, for it may happen to 
be difficult to make J12 exactly equal to J34. 

Regarding the N = 3 chain, the pair density matrix 
in Eq. ^ for the GS shows either a = or d — (for 
Sj = -1/2 or S J = 1/2 respectively). Eq. $Tj$ now 
gives C — 2\z\ = 4|T| for all pairs. We show C12 and C23 
for N = 3 in Fig.[7£a) and (b), respectively. 



B. Longer chains 

For completeness, we consider spin chains extending 
the ones in Fig. [1] following the scheme of Fig. EJa), 
where spin 1 remains at the "surface" . We have veri- 
fied that for longer (A > 4) even-numbered chains, the 
pairs 1 — 2 and 2 — 3 are still exact representations of 
WS. Other pairs may or may not exhibit an entangle- 
ment/disentanglement transition, only with a different 
functional form for C(T) due to the interactions of spin 4 
with additional spins (as compared with the A — 4 chain 
where 1 and 4 are equivalent by symmetry). 

For odd-numbered chains, C12 and C23 for N = 3,5,7 
and 9 are presented in Figs. [T^b) andJTlJc). It is clear 
that the results are qualitatively different for the two 
selected spin pairs. The behavior of C12 is remnant of 
the linear dependence on |T| exhibited by the N = 3 
chains [Fig. E^b)] while C23 for odd N > 5 presents a 
entanglement / disentanglement transition at a value of T 
that approaches the one for a WS (—1/12) for increas- 
ing N, as shown in Fig. [7fc). We infer that zero total 
spin, S J = 0, is not a requirement to obtain entangle- 
ment/disentanglement transitions, even though the in- 
volved states may not be mapped into the standard WS. 



VI. DISCUSSIONS AND CONCLUSIONS 

Our study of quantum and classical pair correlations 
in the GS of open linear AF spin chains devotes partic- 
ular attention to short chains, namely N = 3 and A = 4 
chains. The qualitatively distinct behavior of even and 



odd— A chains, already well established^ is further il- 
lustrated by correlation properties. 

For a class of states in even-numbered chains (here 
fully exemplified by A = 4), we identify two-spin den- 
sity matrices as exact representations of WS, and we 
present the relationship between the WS weight param- 
eter and the spin-spin correlation function of the spin 
pair. In our A — 4 example, we considered a chain with 
a tunable central bond J\ while the others are fixed to 
Jo- We remark that if the "fixed" bonds are different, 
i.e., if the bonds sequence is Jq — J\ — J2, the entan- 
glement/disentanglement transition takes place at a dif- 
ferent value of Ji, which is now a function of J 2- We 
verified that even in this more general case exact WS 
representations are obtained. These results point to a 
possible experimental WS realization in the framework 
of solid state systems, namely the GS of an AFM spin 
chain. This is particularly interesting because so far WS 
realizations were mainly proposed within quantum optics 
frameworks i^ir— We note that our result was obtained 
for zero temperature. It has been shown that the WS 
mixing coefficient may be associated with the tempera- 
ture of a one-dimensional Hcisenbcrg spin models 

Finally we suggest possible physical systems for exper- 
imental verification of WS behavior in spin chains. The 
candidates include: (i) donor-bound electron spins in a 
donors array precisely positioned in Si. Such arrays with 
up to four P donors are within fabrication capabilities, 
as discussed in Ref. [3^ (ii) Quantum dots arrays with 
single electron occupation each. Similar four-dots arrays 
fabrication were also recently reported^; (iii) magnetic 
atoms nanochains4£. In systems (i) and (ii), control over 
Ji is possible through nano-electrodes aligned to the re- 
spective bond. Observation of WS in simple condensed 
matter systems as proposed here seems to be accessible 
with current nano-processing capabilities. 
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